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Many-body localization in an XY model with a long-range interaction is investigated. We show 
that in the regime of a high strength of disordering compared to the interaction an off-resonant 
flip-flop spin-spin interaction (hopping) generates the effective Ising interactions of spins in the 
third order of perturbation theory in a hopping. The combination of hopping and induced Ising 
interactions for the power law distance dependent hopping V{R) oc R~°‘ always leads to the local¬ 
ization breakdown in a thermodynamic limit of an infinite system at a < 3d/2 where d is a system 
dimension. The delocalization takes place due to the induced Ising interactions U{R) oc R~^°‘ of 
“extended” resonant pairs. This prediction is consistent with the numerical finite size scaling in 
one-dimensional systems. Many-body localization in XY model is more stable with respect to the 
long-range interaction compared to a many-body problem with similar Ising and Heisenberg inter¬ 
actions requiring a > 2d which makes the practical implementations of this model more attractive 
for quantum information applications. The full summary of dimension constraints and localization 
threshold size dependencies for many-body localization in the case of combined Ising and hopping 
interactions is obtained using this and previous work and it is the subject for the future experimental 
verification using cold atomic systems. 

PACS numbers: 73.23.-b 72.70.-|-m 71.55.Jv 73.61.Jc 


I. INTRODUCTION 

Many-body localization-delocalization transition 
serves as a natural generalization of a single particle 
Anderson localization concept U to interacting quantum 
systems at finite temperature^, Q and as a quantum 
mechanical extension of a classical transition between 
deterministic and chaotic behaviors i 0 . Localized 
systems show non-ergodic behavior where each small 
subsystem of it remembers its initial state during an 
infinite time. In the delocalized state the system serves 
as a thermal bath for each part of it (e. g. spin or 
particle) leading them to their thermal equilibrium. 
Many-body localization has been considered in a variety 
of physical systems [§-[Tlj| and it attracts growing 
attention because of its significance in quantum infor¬ 
matics (l3 - [l3| . Indeed, in the regime of localization the 
system remembers its initial state infinitely long, while 
in chaotic state this information is quickly erased by 
irreversible dynamics. Recently it has been suggested 
to investigate many-body localization using cold atomic 
systems [ij, fla - H^ , which can emulate various models 
for localization-delocalization transitions. Particularly 
the model considered in the present paper can be 
realized in diatomic alkali systems [i0 . 

Many-body interaction often leads to the single¬ 
particle localization breakdown because it opens new 
channels for energy or particle transport 0, E ini, m, 
[Tsl - i^ (see however Refs. 0, 0, [2^ where the localiz¬ 
ing effect of quasi-static interaction has been exploited). 
Particularly, the long-range interaction decreasing with 
the distance according to the power law due to dipolar, 
magnetic or elastic forces dramatically enlarges the num¬ 
ber of delocalization pathways because of long-distance 


connections. 

Indeed, even at zero temperature excitations are delo¬ 
calized at arbitrarily strong disorder in an infinite system 
if the interaction decreases with the distance slower than 
(see Refs. [I|,[2|). At a finite temperature and in the 
presence of a long-range interparticle (Ising) interaction 
UijSfSj (or UijTiinj for quasiparticles) decreasing with 
the distance as Uij oc the inevitable delocalization 
is expected at /3 < 2d (see Refs. [lllliEllll and the anal¬ 
ysis of delocalization in Sec. IIV Cl) . This strong dimen¬ 
sion constraint results from the resonant energy exchange 
between the flip-flop transitions of spin pairs caused by 
Ising many-body interactions (see Fig. [TJ. However there 
is no such interaction in the XY model containing only 
hopping terms VijSfSj (below we refer to this term as 
a hopping interaction emphasizing its both hopping and 
interaction nature). Therefore it is unclear how the long- 
range hopping would affect the localization there. 

To address this fundamental question and All the ex¬ 
isting gap in the dimensional constraints obtained only 
in the case of dominating Ising interaction [ll|, [T^ (2^ we 
investigate the effect of the long-range hopping interac¬ 
tion on the many-body localization in a random strongly 
disordered XY model for interacting spins 1/2. We show 
that the effective Ising interaction between spins still ex¬ 
ists and it is generated in the third order of perturbation 
theory in the hopping interaction Vij. This interaction 
decreases with the distance as R~^°‘ and it leads to the 
localization breakdown for a < 3(i/2 in agreement with 
the “extended pair” criterion of Ref. [l0| . 

An XY model is relevant for a variety of phenom¬ 
ena including energy transport in Josephson junction ar¬ 
rays p5j| . exciton transport in quantum gases [1113 
and many other problems of interest. This model is also 
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relevant for quantum informatics including quantum in¬ 
formation processing using quantum dot spins [l^ and 
quantum information transport M- Therefore it is im¬ 
portant to characterize the dynamics in this model par¬ 
ticularly understanding the specific of many-body local¬ 
ization there. 

The paper is organized as following. The model is for¬ 
mulated in Sec. im The effective Ising interaction is de¬ 
rived in Sec. nil and the many-body localization prob¬ 
lem in the presence of natural and induced interactions is 
considered in Sec. IlYl The numerical verification of the 
results using the finite size scaling method is reported in 
Sec. El The conclusions are formulated in Sec. ED 


II. MODEL 


In this paper we investigate a many-body localization 
in an XY model of N interacting spins 1/2 

= + (1) 

i i^j 

occupying d-dimensional hypercube with the density 
n. Fields (jji not correlated in different sites i 
and they are uniformly distributed within the domain 
{—Wl2,Wl2). The breakdown of many-body localiza¬ 
tion is investigated for the hopping interaction decreasing 
with the distance according to the power law 




V 


{ndR,jY 


( 2 ) 


where V estimates the hopping interaction at the average 
distance. 

We assume that the disorder is strong 1/ <C W so 
the only long-range interaction can lead to delocaliza¬ 
tion. The consideration is restricted to the case a > d 
so the single-particle delocalization can be neglected (see 
discussion in Sec. lIVDI for the threshold case a = d). Ac¬ 
cordingly there is no mobility edge for the single particle 
states which are all strongly localized (cf. Ref. 0). We 
consider the infinite temperature limit similarly to the 
previous work [13 because it is more convenient for 
analytical and numerical considerations while its gener¬ 
alization to the finite non-zero temperature is straight¬ 
forward ( 13 . 

In contrast with t he p reviously considered models of 
interacting spins [TD, [la [3 with mixed interactions 

^ V^.SfS- + i ^ U.jStSY (3) 

i i=jtj i^j 

Eq. ([T} lacks the Ising spin-spin interaction Uij. In ad¬ 
dition to the pure XY model we will also discuss the 
general model Eq. (|31) in Sec. IIVBI with the special at¬ 
tention to the regime a < jd not considered yet. To make 
our consideration as general as possible we introduce the 


second independent interaction parameter U = Uni (cf. 
Eq. ([5])) representing the Ising interaction at the average 
distance. The large strength of disorder [/ <C IE is also 
assumed for this interaction. 

The Fermion based version of the Hamiltonian Eq. © 
can be generated replacing spin operators S , 
with Fermi-operators aY a, n = a (see Ref. [SO]) as 

H = Y^ ^ Vija \Oj + ^ XI ■ (4) 

i i^j 

The latter model lacks many-body behavior in the ab¬ 
sence of interparticle interaction {Uij = 0) and it can be 
described in terms of independent Fermions which are 
localized in the case of a large strength of disorder for 
a < d. The XY model Eq. (IT|) behaves identically for 
the one dimensional problem with the nearest neighbor 
interaction , but shows different behavior if there ex¬ 
ists triples of spins i,j,k coupled by non-zero interactions 
Vij, Vjk and Vtk as shown below in Sec. IIIII 


III. INDUCED ISING INTERACTION IN XY 
MODEL 

A. Motivation and zeroth order approximation 


Here we show that the effective Ising interaction ex¬ 
ists in the XY model with the long-range interaction. 
It is generated in the third order of perturbation the¬ 
ory in hopping interactions under conditions of a large 
strength of disorder (see Fig. [T] and derivation below). 
This interaction is relatively weak; yet it can lead to the 
delocalization. 

The main reason for that is the dramatic significance of 
the Isin g in teraction for the many-body energy transport 
PE [13 • Indeed the hopping interaction Vij Sf Sj 
bounds only sites with close random fields (j)i « Yi/ 
while the Ising interaction is capable to induce the en¬ 
ergy transport between pairs of such sites i,i' and j,j' 
bound by the resonant condition (pi — (pit Ri (j)j — (j)j' 
while energies (f>i and (pj can be very far from each other 
(see Fig. [T|). This transport channel can lead to the 
breakdown of many-body localization (see Ref. P and 
Sec. IV). In addition it is fundamentally interesting to 
investigate the appearance of the difference between the 
random XY model Eq. o and its Fermion counterpart 
Eq. (|4]). The third order effect that we study is the first 
non-vanishing contribution to this difference. 

To find the induced Ising interaction the perturbation 
theory is developed for the system Eq. ©. One can sep¬ 
arate the Hamiltonian into two parts including the ran¬ 
dom field term dominating in the case of a large strength 
of disorder 


Ho = J2cP,Sp (5) 

i 
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and considered as a zeroth order approximation and the 
hopping interaction term 

V = Y.V,,StS- (6) 

treated as a perturbation. This perturbation can be sepa¬ 
rated itself into resonant and non-resonant parts depend¬ 
ing on whether the change of random field energy (j)i — (pj 
due to flip-flop transition is large or small compared to 
the hopping amplitudes Vij. Under conditions of strong 
disordered system U <C lU one can introduce an inter¬ 
mediate crossover energy such that U <C </>* <C lU 
and treat the interaction part 

/ 

= E ylJsts-0{\cp^ -p,\- p.) (7) 


is created using the transformation 


Heff = exp Hexp , 


(9) 


The matrix in exponent S is chosen to satisfy the condi¬ 
tion [ 3 ^ 


[S,Ho] = -tff. (10) 

In the case under consideration of a strongly disordered 
system one can treat S' as a perturbation and restrict its 
consideration to the few lowest order contributions which 
can be written as (cf. Ref. [s^) 


as an off-resonant perturbation that cannot lead to real 
transitions. Here 0{x) = 1 if a; > 0 or 0 otherwise. This 
part of interaction will be used below to generate the ef¬ 
fective Ising interaction between spins. Since we choose 
4>^ <^W the vast majority of interactions belong to the 
off-resonant perturbation Eq. ©• The remaining reso¬ 
nant interaction 

/ 

= E vl^sts-e{-\p. -cp,\ + P.) (8) 

should be treated as it is. 

In the next step we perform the approximate unitary 
transformation of the Hamiltonian targeted to remove 
off-resonant interactions, which is quite similar to the 
one used in Ref. [3^ ]. This transformation suggests 
the following procedure. The new effective Hamiltonian 


i7e// = i7o + + 17(2) + U(3); 

1 iV 1 


U(2) = i 
2 1 


j/(3) ^ i 
2 1 


S' 


S,Vr 


S,Voff 
3 L 


S,U 


> Voff 


( 11 ) 


We need the third order contribution because only in this 
order the effective Ising interaction of spins is generated. 
The analysis of induced interactions is performed below 
in Sec. IHIBI 

B. Calculation of induced interactions 


1. Second order contributions 


Using Eq. (fTTI) one can expressed the interaction in¬ 
duced in the second order in the off-resonant hopping as 


J 


u(2) = + i/U) ^ . 




-to- 




-,u,s+s- 




U(2) = 


y(2) _ 

U// “ 


-Y^sts-eip. - \4>, - -cPk\- cp.)e{\cp, -M-cp.)- 

-‘^Y^tSj9(.(p* - \(t>i - (t>j\)Y ^ ~ 1 '^* “ (t>k\)9{\(pj - (t>k\- (p*)\ 

i^j k 


i^3 

-2YstS-^{\<P^-h\-^*)Y 

k 


(02 0fc)(0j 0A:) 

(02 0fc)(0j 0fc) 


^( 0 * 10 ^ 0 fc|)^(l 0 i 0 fc I 0 *)- 


( 12 ) 


The results include the corrections to the random field resonant interactions {V^'^ and Voff^ respectively). The 
vY and the second order corrections to resonant and off 
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correction to the random field is much smaller than that 
random held itself. The off-resonant part of interaction 
can be removed applying one more unitary transforma¬ 
tion similarly to Eq. m and it can be used to generate 
the Ising interaction only in the fourth order in the hop¬ 
ping. The resonant interaction represents a weak correc¬ 
tion to the original resonant interaction Eq. m due to 
the large strength of disordering compared to the hop¬ 
ping. Thus all these interactions can be approximately 
neglected. 

2. Induced Ising interaction 

The hrst non-zero contribution to the Ising interaction 
appears in the third order in the hopping interaction from 


three spin loops (see spins i, j and k in Fig. [1] cf. Ref. 
[H). We consider only Ising off-resonant contribution 
to the interaction as the signihcant third order correc¬ 
tion (resonant contribution is smaller as discussed below 
and off-resonant non-Ising contribution can be removed 
by the additional unitary transformation). This contri¬ 
bution is generated by the last term in the definition of 
the third order interaction Eq. (ITT]) and it can be 
expressed as 


- 


'^kS^S- 

yjkS^s^. 

+ ^ V 

'VjkS^S- 

]4i5/'5j c+e-1 

S- 

1 

_1 


3 “ 

J i,j,k 

1 

1 

[Pk-P. 


(13) 


where the notation means that the summation is over 
only off resonant pairs of spins. One can evaluate the 
commutators in Eq. (USD as 

[S+S-, [5+5,-,5+5-]] = 25|(5,^ - 5f). (14) 

Then after the straightforward calculations we come up 
with the following result for the generated Ising interac¬ 
tion of off resonant spins 1/2 within the x — y model 


faster or like R (see Ref. [13), so the sign variable 
sums over spins k in Eq. (US possess the Levy statistics. 
This means that they are determined by shortest possible 
distance of order of the average distance between spins. 

Typical random field differences for few involved spins 
are generally given by the energy disorder range (pi — (pj ^ 
4>i~4>k <Pj — <Pk W. Consequently, the Ising spin-spin 
interaction U-pSpSj Eq. (ITSl) can be estimated as 










{(pi <Pk){<Pj <pk) 


(15) 


There is no such correction in a one-dimensional XY 
model with the nearest neighbor interactions because it 
requires spin pairs i and j, j and k and z, A: to be near¬ 
est neighbors simultaneously which is not possible in a 
one-dimension. In other words the product VijVikVjk is 
always equal zero in that model. Therefore it was neces¬ 
sary to introduce the nearest neighbor Ising interaction 
in Ref. to obtain a many-body behavior. 

The result Eq. (iTSl) is valid only if energy differences in 
denominators are not very small \(pi — pk\ > P*- However 
the opposite situations occur rarely because of the large 
random potential strength W compared to the spin-spin 
interactions. Indeed, the dominating contribution to the 
Ising spin-spin interaction of spins i and j comes from 
the “assisting” spins k located in the direct vicinity of 
either spin i or spin j. This is because the interaction 
is weak {Vtk^ Vjt <C W) and decreases with the distance 




IT2 


(16) 


(remember that V stands for the nearest neighbor inter¬ 
action Eq. ©)■ Consequently the induced Ising interac¬ 
tion at the average distance can be estimated as 


0 ( 3 ) 


W^- 


(17) 


This result will be used to investigate many-body local¬ 
ization in Sec. m occurring due to the combination of 
resonant hopping and off-resonant Ising interactions. Be¬ 
low in Sec. IIII Cl we briefly discuss the induced interac¬ 
tion in other models of interest. 


C. Other models 

Consider the induced interaction for the problem of 
Fermi particles Eq. (|T]). In the case of only non-zero hop¬ 
ping {Uij = 0) there is no induced interactions containing 
the products of two or more operators rii = afui. For 
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instance the third order correction to the energy analo¬ 
gous to Eq. (fT^ does not induce interaction because the 
commutators similar to those in Eq. (TT^ lead only to 
the corrections to a random potential. Indeed one has 

[a+afc,a+aj]] = rii - Uj, ( 18 ) 

and the corresponding correction to random potentials 
can be expressed as 


(a > d) and small pairs (/3 > d) constraints describe the 
stability of the localized state with respect to the single 
particle delocalization 

The third scenario of “extended pairs” (see Fig. [1]) is 
related to large resonant pairs formed by XY interac¬ 
tions. For each resonant pair of spins i, i' [SI -I- SI, = 0) 
the resonant condition [38l| 

( 21 ) 



/ 


E 


Yij Vik Vjk 

[4^i 4^k')[4*j 


(19) 


This is the consequence of the single particle nature of 
the Fermion problem Eq. with Uij =0. 

For the XY model with spin greater than 1/2 the spin- 
spin Ising interaction is induced already in the second 
order of perturbation theory Eq. This interaction 

takes the form 



ij 


2V^S!S^{S^ - S^) 

[fpi — (j)j ) 


( 20 ) 


Consequently it depends on the distance as as for 

the problem with spin 1/2 Eq. (fTBl) . However it is less 
sensitive to disorder. 



FIG. 1: Many body interaction of two resonant pairs of spins 
[i,i') and [j,j') assisted by neighboring spin k. Wave-lines 
indicate significant spin-spin hopping interactions V. 


IV. EFFECT OF INTERACTIONS ON 
MANY-BODY LOCALIZATION IN XY MODEL 

A. Delocalization in XY model due to long-range 
interaction 

Consider the joint effect of induced Ising Uij and orig¬ 
inal hopping Vij interactions on the many-body localiza¬ 
tion using the previously developed theory [ij, [IE Ull in 
the case of the hopping and Ising interactions decreasing 
with the distance according to the power law character¬ 
ized by exponents a and /3 = 2a, respectively and by 
interaction constants given by Eqs. @ and (113- 

Several delocalization scenarios have been considered 
in Ref. [lE| and the critical dimensions (or critical power 
law interaction exponents a and fd) for many-body local¬ 
ization have been suggested (we consider only anisotropic 
interaction case). First two scenarios including hopping 


is satisfied (in general the perturbation theory corrections 
to random fields (j) like Eq. (jT^ comparable to the in¬ 
teraction Vij should be included into field definitions as 
discussed in Refs. [ll|) [2^1. Then two states 5f = 1/2, 
Sf, = —1/2 and Sf = —1/2, Sf, = 1/2 are represented 
nearly equally in the system eigenstates. 

The Ising interaction of resonant pairs [i, i') and [j,j') 
can lead to the energy hopping between them at small en¬ 
ergy differences [ipi — 4>ii) — [4>j — (j)j ') compared to their 
Ising interaction, while the field difference between two 
pairs (pi — pj can be as large as maximum disorder energy 
W (in this case the estimate Eq. (US is well justified). 
This additional channel of the energy transport leads to 
the new critical dimension constraint [lE | d < In 

the case of XY model one can use the induced interaction 
Eq. dm) characterized by the exponent /3 = 2a. Substi¬ 
tuting /3 = 2a into the extended pair constraint d < 
we obtain the new critical dimension constraint for the 
many-body localization in XY model (see also Table |T]) 


id 

y 


< a. 


( 22 ) 


This is the new constraint compared to the previous 
studies [ll|, imii restricted to the regime of a signifi¬ 
cant Ising interaction (/3 < a), where the fourth scenario 
of “Iterated pairs” determines the critical dimension con¬ 
straint 2d < jd [a (anisotropic interaction is assumed). 
In the case of dominating hopping interaction a < jd the 
extended pair scenario gives the stronger restriction of 
the critical dimension. 

Below we extend the previous consideration to the gen¬ 
eral problem Eq. in the case a < (d and obtain the 
dimension constraint for many-body localization (Sec. 
IIVB|) for arbitrarily relationship of power law interac¬ 
tion exponents a and jd. We also estimate the minimum 
system size where the localization is still possible at di¬ 
mension exceeding the critical one (Sec. IIV C|) . The lat¬ 
ter study is important for real systems having finite size 
where the proposed theory can serve as a background 
for the analysis of the future experimental data. These 
systems can be realized using the chains of cold atoms 
dEilE)- The results are summarized in Tables HI and Ull 


B. Dimension constraint for dominating hopping 
interaction [a < P) 

Here we discuss the dimension constraint for the gen¬ 
eral problem Eq. © with a < jd not addressed previ- 
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TABLE II: Critical strength of disorder and number of particles determining localization threshold for violated dimension 
constraints Table [I] 


Model 


Specific case 


= nRi 


W, 


d < /3 < a, 2d 


a < 


dj3 

/3-d 


^ d(c, + j3)-a/3 d(a + li)-afi 


2d 

2d-/3 


2d 

2d-0 


V -max 


(f) 


d(c. + l3)-c.0 

N - ^ -. 


(«) 


Q 2d — 

Ana 


d0 
0 — d 


< a 




V 




2d-/3 

N^2T- 


d < a = P 


P < 2d 


d'^ d 

Y \ 0(2d-0) f W\ ‘2d-0 


M 


V'ISI’a'^ 


M 


d < a < P, ^ 


P<-^ 

' a — d 


^ L j d(Q + /3)-c</3 


dp 

d(a + P)-aP 


•(v) 


2d(d+a) 

a{3d-2a) 


- ‘ld(d-\-a) 

a:(3d-2a:) 
V 


V -max 


(« 


d(a + P)-aP a(3d-2a} 

jV-33- N 2d(d+c) 


< 


a(3d-2a) 
^ 2d(d + Q;) 


^ <a < P 


P <-^ 

' oc — d 


£\ d{a + p-)-aP / d{a + P)-aP 


, -3 d(c + 3)- 

V" ^ N 33- 




TABLE I: Dimension constraints for many-body localization 
for A —>■ oo 


Model 

a > P 

a < P <2a 

2a <P 


P 

a.0 

2o 


2 

ot + 8 

3 

dcF 

p 

0.0 

a0 

2 

__ 

0+3 


ously. The induced Ising interaction Eq. (HH) can be¬ 
come more significant than the initial Ising interaction 
if it decreases with the distance slower than the original 
Ising interaction i. e. at 2a < p. Therefore consider¬ 
ing the dimension constraint one should characterize the 
Ising interaction by the minimum power law exponent 

P^, = min(/3, 2a), (23) 

and use this exponent in the dimension constraint for 
extended pairs d < [l^ as it is given in Table HI 

This result does not depend on the spin since the induced 


Ising interaction Eq. (j20p for large spins shows the same 
distance dependence as for the spin 1/2 Eq. ([T6)l . 

Eor the Fermi-particles counterpart model with mixed 
interactions (a < /3) Eq. (|4]) there is no induced many- 
body interaction so the extended pair criterion should be 
always applicable (see Table |l|. 

In the case of violated dimension constraint the delo¬ 
calization should take place at arbitrarily large strength 
of disorder beginning from the sufficiently large system 
size. We estimate the maximum number of spins where 
many-body localization can still take place in Sec. IIV Cl 
below. Although this derivation essentially repeats the 
arguments of previous work [HI . [13 there is some 
qualitative difference of the case a < P compared to the 
previously considered case a = /3 so it can be useful for 
better understanding of many-body localization break¬ 
down in this regime and for the analysis of the future 
experimental data in cold atomic systems. 


C. Critical size and disorder strength at large 
dimension d > dc 


We begin with the analysis of delocalization in the gen¬ 
eral case a < P including XY model. Assume that the 
dimension constraint for extended pairs d < is not 

satisfied so the delocalization should take place at ar¬ 
bitrarily large strength of disorder at sufficiently large 
system size. Our targets are the dependence of criti¬ 
cal disorder Wc on the system size L expressed through 
the number of spins N = nL‘^ and the dependence of 
the number of spins Nc needed for the delocalization to 
occur at the given strength of disorder W. Our consid¬ 
eration follows the previous analysis of similar problems 

for a > /3 [Ullil. 


The dynamics of strongly disordered system is primar¬ 
ily associated with resonant pairs Eq. m where the 
random field energy can be compensated by the hop¬ 
ping interaction. The resonant pairs can be divided 
into sub classes of typical size R and typical energy 
V{R) = —r-—■ The density of resonant pairs having 

{nSR)<=‘ 

typical size R is given by [22l| 


n{R) 


V 1 
n -5-. 

W (ndR)^-d 


(24) 


Ising interaction characterized by the interaction 
strength at the average distance U creates the dynamic 
interaction of resonant pairs of the typical size R lead¬ 
ing to the simultaneous flip flop transition in both reso¬ 
nant pairs accompanied by the energy hop between them. 
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Since both pairs are resonant the amplitude of this simul¬ 
taneous transition is of order of their bare Ising interac¬ 
tion which can be estimated as [1^ 


U{R) = Un{R)^U 



1 




3 


(25) 


One should notice that in the case of extended pairs and 
a < 2d the average distance between pairs n(R)~d be¬ 
comes smaller than their size at sufficiently large size R 
(see Ref. [l^, Eq. (l24l) l for two resonant pairs. There¬ 
fore the coupling of only the closest spins is important 
as shown in Fig. [TJ Then the isotropic or anisotropic 
character of Ising interaction does not affect the estimate 
of Ising interaction between pairs Eq. (ESI) and conse¬ 
quently the delocalization transition in contrast with the 
case a > /3 of dominating Ising interaction considered in 
Ref. 0. 

The interaction of resonant pairs Eq. E3 represents 
their dynamic coupling that can lead to delocalization if 
it exceeds the typical energy V{R) of the subsystem of 
pairs. Considering power law distance dependencies of 
U (R) and V (R) one can see that if the criterion d < 
is not satisfied then the interaction of pairs of the size 
R {U{R)) decreases with the distance slower than the 
typical energy of those pairs {V{R)). Therefor inevitable 
delocalization is expected at sufficiently large system size 
where the pair coupling exceeds their typical energy. The 
size R* where two energies are equal to each other esti¬ 
mates the maximum system size where many-body local¬ 
ization is still possible at the given strength of disorder 
W. Solving the equation C/(i?*) = V(R^) we estimate 
this size as 


R* =n 


d(c + ^)-cf) 


(26) 


In the case of XY model with the spin 1/2 the Ising 
interaction strength U depends itself on the strength of 
disorder Eq. (fTBl) . In this case the dependence of the 
maximum system size on disorder Eq. (j26p should be 
modified as (remember that /3 = 2a in that case) 


R*xy 


_ 2 ^ 
= n d. 



2(d+oc) 

a{3d-2ac) 


(27) 


Consider the general model with mixed interaction Eq. 
(El for a < 13. If both constraints for mixed interactions 
(a/3/(a + (3) > d) and for XY model (a < 3d/2) are vio¬ 
lated then one should choose the strongest interaction out 
of two leading to the lowest size constraint determined by 
the minimum of two estimates in Eqs. Ei) and EH)- If 
only one of two conditions is satisfied then one should use 
the corresponding criterion. The estimates of the critical 
number of spins N = nR'^ in all these cases is given in 
Table m 

Case a = P considered in Ref. can be described by 
Eq. (05)) setting a = p. In the remaining case oi a < P 


the iterated pair scenario has been suggested [l^ (see also 
[llj). In that scenario the delocalization is determined by 
interacting nearest neighbor resonant pairs which replace 
spins in the original model. For anisotropic interaction 
|16l | the problem can be reformulated with the modified 
parameters Wp = V,np = nV/W, Up = Vp = U{VjW)^ 
and Up = Pp = p. Using these parameters in Eq. Ei) 
one can estimate the critical size assuming P < 2d as 


R*p 




(28) 


If both extended pair and iterated pair criteria lead to 
the delocalization {a > P, P < 2d, aP/{a + P) < d) 
then one should choose the minimum size determined by 
either Eq. (l26l) or Eq. (l28l) corresponding to the most 
efficient delocalization. 

In TableUwe show the estimates of the maximum num¬ 
ber of spins {N^, = nR'p) where many-body localization 
is still possible for all possible exponents a and P and 
violated dimension constraints. The critical strength of 
disorder Wc for the given number of spins can be found 
resolving Eqs. (1^ . (I?7l) . (E51) with respect to the disor¬ 
der W at fixed system size R or number of spins X. 
The dependencies Wc{N) are also shown in Table [ill 
The results can be easily generalized for the interact¬ 
ing Fermions simply removing the parts related to the 
induced Ising interaction which does not exist in that 
case. 


D. Threshold regimes d = 

The threshold regimes where interactions behavior cor¬ 
responds to the border line of the dimension constraints 
need the special attention. There are two distinguish¬ 
able regimes including the cases a = d or P = d, where a 
single particle delocalization becomes important at very 
large system size for arbitrarily large strength of dis¬ 
order [i^, [23 and the cases P = 2d, a = 3d/2 and 
aPj{a + P) = d corresponding to the border line for 
many-body localization breakdown. The first regime is 
particularly significant because of the common appear¬ 
ance of R~‘^ interactions due to dipolar or elastic forces 

We believe that the single-particle delocalization can 
be ignored for a = d (and similarly for P = d) in the 
regime of strong disorder U <C lU because the delocal¬ 
ization becomes significant only if the logarithmic de- 
localization parameter In (IV) approaches unity [^ . 
This can happen only at exponentially large system size 
remarkably exceeding the critical sizes from Table HIl for 
many-body localization. Therefore we do not expect any 
related changes in our estimates for critical size and dis¬ 
order behaviors except for may be additional weak log¬ 
arithmic parametric dependencies which are beyond the 
scope of our qualitative analysis. 







The number of resonant interactions grows logarithmi¬ 
cally with the system size at the many-body localization 
thresholds (/3 = 2d, a = M/2 or al3/{a + /]) = d). This 
behavior is similar to the single particle localization prob¬ 
lem with strong disorder and R~‘^ hopping interaction 
M- In this regime the excitations are nearly localized: 
their typical localization radius is much smaller than the 
system size, but it approaches infinity in the infinite sys¬ 
tem limit (see however Ref. [s^ where delocalization 
can take place under certain conditions at exponentially 
large but finite radius). In our case we similarly expect 
that the system remains localized until the system size 
gets exponentially large ln(n'si?) ^ WjV while we are 
not able to describe the system behavior at larger sizes. 
Perhaps the experimental investigation of the implemen¬ 
tations of this model in cold atomic gases as suggested in 
Refs, [l^ [l^ will help to resolve this difficult issue. 

Below in Sec. El the predictions of theory for the XY 
model are verified for one-dimensional model using the 
numerical finite size scaling of the ergodicity parameter 

[H. 


V. FINITE SIZE SCALING FOR XY MODEL 


Consider a random one dimensional XY model. 
In the numerical study we describe the localization- 
delocalization transition using the ergodicity parameter 
[12, . In the thermodynamic limit of infinite system 

this parameter is defined as the local spin-spin correlation 
function 


Q = lim 

t—^OO 


< m\S/{t)S/{0)\m > 6{Em) 


(29) 


taken in the infinite time limit [4nl | and averaged over the 
system eigenstates having zero energy. The latter choices 
correspond to the infinite temperature thermodynamic 
limit [ 23 , . In that limit the average ergodicity param¬ 

eter should approach zero in the delocalized state where 
correlations are subject to decay, while in the localized 
state it should be finite (unity in an infinitely strong dis¬ 
order limit). Our method is perhaps not as efficient in the 
definition of the localization transition point as recently 
developed analysis of entanglement entropy [12, lil| ; yet 
it permits the easy determination of the scaling of the 
localization transition with the system size that can be 
compared to the theory predictions in Table [ITl 

The numerical study is performed for the Hamiltonian 
Eq. O describing the periodic one dimensional chain of 
N spins separated by the unit distance using the periodic 
interaction 


= ± 


max(\i 


1 

j\,N - \i-j\)°‘ 


(30) 


with random, uncorrelated signs of all interactions. This 
model is similar to the models studied in Ref. [il,[13 
with the only difference that the Ising spin-spin interac¬ 
tion is set equal to zero here. The total spin projection 


to the z axis is conserved. We restrict the numerical con¬ 
sideration to the states with zero total spin projection 
which represents well the thermodynamic limit. 

In the numerical study the ergodicity parameter Eq. 
(l29l) is defined as a configuration averaged spin-spin cor¬ 
relation function at infinite time (40j | over the narrow 
band of eigenstates a of the Hamiltonian Eq. (ED 
with zero total spin and energy belonging to the do¬ 
main (—6 < Ea < d) 

Q = (31) 

^ a 

The domain size <5 = Q.Q4W^/N is determined requiring 
the many-body density of states 


9{E) 


( 2iNW'^ ) 

^/2AttNW‘^ 


(32) 


to change at the scale 5 by 1%. The function g{E) has 
been estimated using the random potential part of the 
system Hamiltonian Eq. ([SD employing the law of large 
numbers {N ^ 1) and large strength of disordering V <C 

W. 


0.75 


O 0.5 


0.25 



150 


fi/4(14)W/f^^4(N) 


FIG. 2: Original (inset) and rescaled dependencies of ergod¬ 
icity parameter Q on disordering W for 0 = 1. 


The parameter Q has been averaged over few thou¬ 
sands realizations of the system Hamiltonian to have the 
relative error of 1% (5% for N = 16). The XY mod¬ 
els with interactions characterized by power law expo¬ 
nents a = 1,1.25,1.5,1.75, 2 taken around the predicted 
threshold a = 3/2 (see Table HD have been considered 
for a strength of disorder W ranging between W = 1 
and W = 150 (we set V = n = 1). In all systems the 
localization transition can be seen with increasing disor¬ 
der; i. e. Q approaches 0 at small W and tends to 1 in 
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the opposite limit of hK —>■ oo (see Figs. [2l [3l SI [5l [6]) . 
However these transitions behave differently for different 
exponents a and numbers of spins N. 

The results for the dependence of ergodicity parameter 
on the strength of disorder are shown in the insets of Figs. 
SlSlSlfor a < 3/2 and in Figs. [5]and[n]for a = 1.75 and 
2, respectively. For a < 3/2 the transition shifts towards 
large disorder W with increasing the number of spins, 
while for a > 3/2 the transition is almost insensitive to 
the system size. The situation is inconclusive for a = 
3/2 where larger system sizes should be probed. The 
result of the visual inspection agrees with the qualitative 
expectation for the threshold power law exponent ac = 
3/2 separating the regimes where the delocalization is 
inevitable in the large N limit (a < 3/2) and where it is 
possible. 



FIG. 3: Original (inset) and rescaled dependencies of ergod¬ 
icity parameter Q on disorder W for a = 1.25. 

For quantitative characterization of the critical 
strength of disorder Wc dependence on the number of 
spins N we use the scaling function (2^ 

This function replaces the continuous power law depen¬ 
dence iV“ predicted by the qualitative consideration (Ta¬ 
ble ini) with the discrete sum of resonance probabilities 
which accounts better for finite size effects. For a > 0 
this function has the asymptotic behavior fa{N) oc N°‘. 

In Figs, m El m we show the modified plot of the data 
for a < 3/2 with disorder strength rescaled using the 
function Eq. (1331) with the theoretically predicted expo¬ 
nents a = {a — 3d/2)/(l -|- a) (see Eq. (P71) and Table 
El). This rescaling leads to a reasonable match between 
the graphs for different numbers of spins N indicating 



FIG. 4: Original (inset) and rescaled dependencies of ergod¬ 
icity parameter Q on disorder W for a = 1.5. 



FIG. 5: Dependence of ergodicity parameter Q on disorder 
W for a = 1.75. 


a good agreement of theoretically predicted and numeri¬ 
cally found scaling of the transition point with the system 
size. 

In addition to the visual inspection we performed 
quantitative estimates of rescaling parameters for the 
data sets Qn{W) with the same exponent a and differ¬ 
ent N's using the optimization procedure [1^. For each 
exponent a we determined the set of optimum rescaling 
parameters Ca{N) corresponding to the minimum of the 
squared deviation ~ QN{ca{N)Wi)Y with 

Wi changing with the step of 1 from 1 to 100. The re¬ 
sults for different rescaling parameters c{a, N) are shown 
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FIG. 6: Dependence of ergodicity parameter Q on disorder 
W for a = 2. 



FIG. 7: Scaling factors c vs. a number of spins N for different 
interactions (a) and the fit of c{N) dependencies using func¬ 
tions fa{N) with optimum fitting parameters a estimating 
scaling exponents for localization transition size dependence. 
Inset shows the dependence of exponents a on the interaction 
power law exponents a. The solid (red) line indicates the 
threshold case a = 0. 


in Fig. [7] by symbols defined within the graph. 

To characterize the change of the rescaling factors and 
consequently the critical strength of disorder Wc with 
the number of spins we fitted each data set by the func¬ 
tion fa{N) Eq. (l33l) choosing the parameter a to attain 
the best agreement of the model with the numerical re¬ 
sults. The solid lines in Fig. [7] show these optimum fits 
and the exponents a are shown in the inset and indi¬ 


cated in the legend. The scaling function fo{N) serves 
as the crossover between the regimes of finite (a < 0) and 
infinite (a > 0) localization thresholds in the thermody¬ 
namic limit fV —>■ oo. In agreement with the theoretical 
expectations of the crossover at a = 3/2 we found a > 0 
for a = 1,1.25, a < 0 for a = 1.75, 2, while situation is 
not clear at the threshold a = 1.5 as shown in the in¬ 
set to Fig. [71 The estimated exponents a are, however, 
somewhat larger than the theory predictions Eq. (Ell) 
which can be due to finite size effects. The logarithmic 
increase of the number of resonant interactions with in¬ 
creasing the number of spins for a = 1 (see Ref. [2^ and 
Sec. IIVD|) can lead to the observed deviations. 


VI. CONCLUSION 

We investigated many-body localization in the XY 
model with the long-range interaction V(R) oc R~°‘ in 
a strongly disordered regime. In this regime the Ising in¬ 
teraction decreasing with the distance as R~^°‘ appears 
in the third order of perturbation theory with respect 
to the hopping interaction. Considering the combined 
effect of this induced interaction and the original hop¬ 
ping interaction we found the novel dimension constraint 
a > 3d/2 required to attain the many body localization 
in the thermodynamic limit of infinite system. 

Using this result we suggested the dimension con¬ 
straints for the general problem of many-body localiza¬ 
tion in the presence of the long-range hopping and Ising 
interactions with arbitrarily power law distance depen¬ 
dencies (see Table HI. In the case when the dimension 
constraint is violated the maximum system size where 
the localization is still possible is predicted (see Table 
HB. Also we predicted the critical disorder dependence 
on the systems size. These results can be used to inter¬ 
pret the recently proposed experiments in cold atomic 
systems implementing the disordered spin systems with 
the long-range interactions [H, HBl • 

It turns out that many-body localization in XY model 
is more stable with respect to the long-range interaction 
then in the model of s pin s having both Ising and XY 
long-range interactions |22| . For instance many-body lo¬ 
calization can be attained in a thermodynamic limit in a 
three dimensional XY model with a quadrupole interac¬ 
tion (a = 5) while it is not possible in the spin system 
with a quadrupole Ising and XY interactions. Thus spin 
systems with only XY interaction can be more attractive 
for quantum informatics. 

It is interesting to notice that the effective long-range 
interaction of spins can b e g enerated in periodically 
driven interacting systems [Ifl l42l - l43 | as a part of ef¬ 
fective Hamiltonian as suggested in Ref. [i^. Therefore 
the scaling behaviors similar to the ones studied in the 
present work can be seen in these systems as well (e. g. 
the scaling of the decay time with drive parameters |43|). 
The detail analysis of the relationship between two prob¬ 
lems is beyond the scope of the present work and will be 
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